Reconstructing directed networks for better synchronization 
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In this paper, we studied the strategies to enhance synchronization on directed networks by 
manipulating a fixed number of links. We proposed a centrality-based reconstructing (CBR) method, 
where the node centrality is measured by the well-known PageRank algorithm. Extensive numerical 
simulation on many modeled networks demonstrated that the CBR method is more effective in 
facilitating synchronization than the degree-based reconstructing method and random reconstructing 
method for adding or removing links. The reason is that CBR method can effectively narrow the 
incoming degree distribution and reinforce the hierarchical structure of the network. Furthermore, 
we apply the CBR method to links rewiring procedure where at each step one link is removed 
and one new link is added. The CBR method helps to decide which links should be removed or 
added. After several steps, the resulted networks are very close to the optimal structure from the 
evolutionary optimization algorithm. The numerical simulations on the Kuramoto model further 
demonstrate that our method has advantage in shortening the convergence time to synchronization 
on directed networks. 
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I. INTRODUCTION 



Synchronization is an important dynamical process in 
many systems. Understanding and controlling this col- 
lective dynamics is of both theoretical and practical sig- 
nificance [l|, [|| . Many methods have been proposed to 
enhance the network synchronizability, including the re- 
distribution of the coupling strengths [3j-|8| , the modifica- 
tion of the network structure [9l-fl3j , the flipping of link 
directionality p^4l7| . and so on. Each group of meth- 
ods has its specific range of applications since different 
systems are under different technical constrains. 

As real systems are frequently manipulated by grow- 
ing out new connections and eliminating redundant con- 
nections [jq . we here focus on the link reconstructing 
method to enhance synchronizability. The word "re- 
constructing" stands for manipulations including adding, 
removing and rewiring links, and in each manipulation 
only one link is created or changed. The reconstructing 
method for undirected networks has already been studied 
by usingthe information embedded in Laplacian eigen- 
vectors Q- However, the spectral method cannot be di- 
rectly applied to directed networks, since the complex 
value emerges in eigenvectors when the Laplacian matrix 
is asymmetry. Up to now, how reconstructing methods 
affect the synchronization on directed networks lacks of 
systematic consideration. 

Previous works show that there are two important fac- 
tors that mainly affect the synchronization on directed 
networks: the in-degree distribution Q and the hierar- 
chical structure |19l - [23 ]. Generally speaking, the more 
homogeneous the in-degree distribution in a directed net- 
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work is, the stronger the network synchronizability is. 
The synchronizability in directed networks can be ap- 
proximately expressed by R w k^/k^f n [3]. A smaller 
R indicates a better synchronizability. Accordingly, to 
enhance synchronizability in directed networks, we can 
either decrease fc^° x or increase k^f n which are respec- 
tively corresponding to removing the links that point to 
the nodes with maximum in-degree or adding links that 
point to the nodes with minimum in-degree. Besides the 
in-degree distribution which plays the dominant role in 
determining the synchronizability of directed networks, 
the hierarchical structure may also have considerable ef- 
fects [23| . If the reconstructing strategy only aims at 
the homogeneous in-degree distribution, an effective hi- 
erarchical structure cannot be formed and the resultant 
synchronizability may fall into local optimum. There- 
fore, how to choose the starting points of the adding 
and removing links is very important which may affect 
the overall hierarchical structure, and thus further af- 
fects the synchronizability. In this paper, we employ the 
PageRank algorithm [24| to characterize the centralities 
of nodes on directed networks and design a centrality- 
based reconstructing (CBR) method to enhance synchro- 
nizability accordingly. Generally, the information of the 
hidden hierarchical structure is embedded in the node 
centrality gradient 25]. The nodes with high PageR- 
ank scores are more likely to be at the high layer of the 
network, while those with low PageRank scores proba- 
bly locate at the low layer. The basic idea of the CBR 
method is that the newly added links should start from 
high-layer nodes (with high PageRank scores), while we 
should choose the links starting form the low-layer nodes 
(with low PageRank scores) to remove. 

Performing the linear stability analysis of synchroniz- 
ability [26l428l| in directed networks, we find that the CBR 
method is more effective in facilitating synchronization 
than the degree-based reconstructing (DBR) method and 
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random reconstructing (RR) method. Furthermore, we 
apply the CBR method to links rewiring procedure. At 
each step we first remove one link from the network and 
then add a new one according to the removing and adding 
strategies of the CBR method. After several steps, we 
find that the resulted networks are very close to the op- 
timal structure from the evolutionary optimization algo- 
rithm 1 151 . The numerical simulations on the Kuramoto 
model [29j further demonstrate that our method has ad- 
vantage in shortening the convergence time to synchro- 
nization on directed networks. 



II. CENTRALITY-BASED RECONSTRUCTING 
METHOD 



Consider a directed and unweighted network with N 
nodes and E links. Denote by A the adjacency matrix 
where the element Aij — 1 if there is a directed link from 
node i to node j, otherwise 0. Specifically, i is the start 
and j is the end. The average degree of the network 
is k = 



the node i's centrality score at time t (t > 1) is given by 



N ■ 



In the synchronization process, each node 
stands for an oscillator and each directed link represents 
the influence from one oscillator to another. The syn- 
chronization in directed networks is actually formed by 
a top-down centralized control mechanism [19| . The os- 
cillators with higher dynamic centrality are supposed to 
be more influential in driving the whole system to syn- 
chronize. In directed networks, there are many different 
ways to measure the node's centrality [24l. [30r 33] . One 
of the prominent group of methods is based on the ran- 
dom walk process, such as PageRank [24[ , HITs [Hf and 
LeaderRank [33[ algorithms. In this paper, we apply the 
well-known PageRank algorithm to quantify the node's 
centrality in directed network. 

PageRank is a famous ranking algorithm which forms 
the basis of Google™ search engine [24{. It has been 
applied to rank scientists in citation networks 13411 and 
detect community structure in directed networks 35] . In 
practice, PageRank assigns a score s, to denote the at- 
tractiveness of the webpage i. Webpage i obtains higher 
score if many other important webpages point to it. From 
the physical perspective, PageRank describes a random 
walk process on directed network, where the score s$ is 
proportional to the frequency of visits to a particular 
node i by a random walker. In PageRank algorithm, 
a parameter c, called return probability, is introduced, 
which represents the probability for a random walker to 
jump to a random node, and 1 — c is the probability for 
the random walker to continue walking through the di- 
rected links. Note that, in the oscillator network the link 
from node i to node j indicates that node i has influ- 
ence on node j, and the nodes who has larger influence 
on other nodes are more important for the network syn- 
chronization. Therefore, in our case the random walker 
should follow the opposite direction of links. In this way, 
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where S a ,b = 1 when a — b and 5 a .b = otherwise. Ini- 
tially, we assign each node one random walker, namely 
Si(0) = 1 for i = 1, 2, • • • , N. The typical value of return 
probability is about 0.15 24j. The final score of each 
node is defined as the steady value after the convergence 

Of Si(t). 

Generally speaking, the hierarchical structure of a di- 
rected network can be embodied by the node's centrality 
measured by Eq. ([!} (25[. A node with higher score is 
inclined to locate at the high layer of the network hier- 
archy. With this in mind, we propose a centrality-based 
reconstructing (CBR) method to enhance the synchro- 
nizability on directed networks by manipulating a fixed 
number of links. Since the network synchronizability can 
be characterized by the indicator R « fcmax/^min : the 
smaller R, the higher synchronizability. Therefore, we 
have two ways to enhance synchronizability: increasing 
A-min by adding links pointing to the nodes with minimum 
in-degree and decreasing fc*" x by removing links point- 
ing to the nodes with maximum in-degree. The starting 
nodes for adding and removing links are selected accord- 
ing to the nodes' centrality scores given by Eq. ([T]). To 
add a link, we choose node i with minimum in-degree 
as the end, and node j with maximum centrality among 
all the nodes having not yet pointed to i as the start. 
To remove a link, we firstly guarantee the end is of maxi- 
mum in-degree, and then among all possible candidates of 
starts, we choose the one with minimum centrality. This 
strategy will help to generate more receptors, namely the 
nodes without any outgoing link, which is a favorable 
structure for synchronization. For details, one can see 
the discussion in Sec. Ill B. Note that, after adding or 
removing a link the centrality score of each node will be 
recalculated. 

We consider two other methods for comparison: the 
random reconstructing (RR) method and the degree- 
based reconstructing (DBR) method. The strategies of 
choosing the ends of these two methods are the same as 
CBR, namely adding new link pointing to the node with 
minimum in-degree and removing link that points to the 
node with maximum in-degree. The difference is how to 
choose the starts. When adding or removing links by the 
RR method, the start of the links are randomly selected. 
For DBR method, when adding a link, the node with 
maximum out-degree and not yet pointing to the end 
will be chosen as the start, while when removing a link, 
the link starting with a node of minimum out-degree is 
removed. In the following section, we will compare CBR 
method with RR and DBR methods on three directed 
network models based on the linear stability analysis of 
synchronizability and the numerical simulation on the 
Kuramoto model. 
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FIG. 1: (Color online) The performance of the RR, DBR and 
CBR methods on synchronizability when adding links to (a) 
directed regular networks (N = 100, k = 1), (b) directed WS 
networks (N = 100, k = 3, p = 0.1) and (c) directed BA 
networks (N = 100, No = 10, fcmax = 6). The results are 
averaged over 100 independent realizations. 



III. RESULTS 

Under the framework of master stability analysis, the 
synchronizability of an undirected network can be quan- 
tified by the eigenvalue ratio of the corresponding Lapla- 
cian matrix of this network, namely R = where 
Xn and A2 are respectively the largest and the smallest 
non-zero eigenvalues of the Laplacian matrix |26-[2c| . In 
directed networks, since the Laplacian matrix, denned as 
Lij = kl n Sij — a,ij, is asymmetric with zero rowsum and 
has complex eigenvalues. In order to achieve the synchro- 
nization condition, every eigenvalue should be entirely 
contained in the region of negative Lyapunov exponent 
for the particular master stability function. If the stabil- 
ity zone is bounded and the imaginary part of complex 
eigenvalue is small enough, the network synchronizabil- 
ity can be approximately measured by the real part of 
eigenvalue ratio R = A^/A^ , where A^ and A2 are re- 
spectively the larg est and the second smallest real parts 
of eigenvalues |l9l I20I ] . Generally speaking, the smaller 
the ratio R, the stronger the synchronizability. Here, we 
employ the indicator R to evaluate the synchronizability 
of directed networks. 

In this paper, we mainly consider three kinds of di- 
rected networks: (i) Directed regular network with iden- 
tical degree k and clockwise links. Here k indicates ei- 
ther the average in-degree or the average out-degree since 
these two values are the same. For example, k = 1 
means that each node has one in-link and one out-link, 
(ii) A variant of Watts-Strogatz (WS) network [H (di- 
rected WS network). Starting from a directed regular 
network, each link will be reconnected with two randomly 
selected nodes with probability p € (0, 1). (iii) A variant 
of Barabasi- Albert (BA) network [37| (directed BA net- 
work). Starting from a directed tree with Nq nodes, a 
new node with m links (m is a random integer between 
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FIG. 2: (Color online) (a) The obtained network by adding 50 
links to a directed regular network with the RR method, (b) 
The obtained network by adding 50 links to a directed regular 
network with the CBR method, (c) The dependence of core- 
distance on the number of the added links in the directed WS 
networks, (d) The dependence of the number of inverted links 
on the number of the added links in the directed BA networks. 
The parameters of the network models are the same as those 
in Fig. [1] The results in (c) and (d) are averaged over 100 
independent realizations. 



1 and fc max ) is added to the network in each step un- 
til the total node number reaches N. Each new added 
link connects to an existing node i with the probability 

Qi = ki ^f' 1 ■ The link direction is set to be from older 
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nodes to younger nodes. 



A. Adding 

The results for adding links in three modeled net- 
works are reported in Fig. 1. Generally speaking, CBR 
performs best among all three methods. Specifically, 
in directed regular networks [see Fig. 1(a)] all these 
three methods can decrease R by adding links. The en- 
hancement of synchronizability results from the decreas- 
ing of the average shortest distance of networks. As we 
know, networks with small average shortest distance are 
likely to have better synchronizability [H, [39J. Besides, 
the DBR method and CBR method outperform the RR 
method especially when many links are added. The rea- 
son is that the DBR method and CBR method are capa- 
ble to generate a network core during the adding process, 
which may involve one (when adding a small number of 
links) or several (when adding many links) nodes. The 
core node has many outgoing links and drives other nodes 
during the synchronizing process. It has been pointed out 
that such a core is favorable for the network synchroniz- 
ability [15]. In contrast, the RR method is unable to 
induce such an effective core, thus its synchronizability 
enhancement is only from the decreasing of the average 
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shortest distance. Two typical examples are shown in 
Fig. 2(a) and Fig. 2(b) which are obtained by adding 50 
links to a directed regular networks with RR and CBR 
methods, respectively. Clearly, there exists a center node 
as the core of the obtained CBR network, yet there is no 
observable central nodes for the obtained RR network. 

Figure [Tfb) shows that the indicator R decreases with 
the increasing number of added links for all these three 
methods and the difference between them is very small 
(-Rcbr = 6.81,i?DBR = 7.02,i? RR = 7.16). The enhancement 
mainly comes from narrowing the in-degree distribution. 
For the directed WS networks, even though the CBR and 
DBR methods can still form a core by adding a few links, 
the positive effect of the core is disrupted by the compli- 
cated structure of the initial network. Therefore, when 
adding a few links, like 50, comparing with RR method 
the synchronizability enhancements by CBR and DBR 
methods are very limited, respectively 4.9% and 2.0%. 
Note that, when we add more links, like 500, the im- 
provements increase to 12.3% and 9.7% for the CBR and 
DBR methods (i? CB R - 1.71,-R DB r = 1.76,i? RR = 1.95), re- 
spectively, since the effect of the initial network will be 
depressed. Additionally, the CBR and DBR methods can 
lead to shorter convergence time than the RR method due 
to the smaller core-distance (see details in Sec. Ill D for 
Kuramoto model). Given a network, the core-distance 
is defined as the average shortest distance from the core 
node to all other nodes in this network. It was pointed 
out that the smaller the core-distance, the quicker the 
whole network converges to synchronization, since the 
core actually drives the whole network to the synchro- 
nized state [l5| . The dependence of the core-distance on 
the number of added links in the directed WS networks 
is shown in Fig. [21(c). Note that, since there is no core in 
RR network, we use the minimal distance from a node to 
all other nodes in the network (such node is considered 
as a core) . Clearly, the values of the core-distance of the 
CBR and DBR methods are much smaller than that of 
RR, and thus the CBR and DBR methods can result in 
shorter convergence time. 

In directed BA models, as shown in Fig.[T](c) it is clear 
that CBR method performs best. The directed BA model 
has obvious hierarchical structure where the nodes with 
higher centrality scores incline to locate at the higher 
layer. In CBR method, the start of the new link is cho- 
sen as the node with maximum centrality score, thus the 
generated core is actually the root node (in the highest 
layer) of the initial directed BA networks. In contrast, 
with DBR method the core nodes usually locate in low 
layers. With starts in low layers many inverse links (i.e., 
links from low layers to high layers) will be formed, lead- 
ing to a large number of loops. The situation is even 
worse with RR method. The dependence of the number 
of inverse links on the number of added links of BA net- 
works is shown in Fig.[2]Jd). We can see that the number 
of inverse links generated by CBR method is fewer than 
DBR method and RR method. The large number of loops 
are not beneficial to synchronization [lg . |40| , and thus in 
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FIG. 3: (Color online) The performance of the RR, DBR and 
CBR methods on synchronizability when removing links from 
(a) directed regular networks (N = 100, k = 3), (b) directed 
WS networks (N = 100, k = 3, p = 0.1) and (c) directed BA 
networks (TV = 100, No = 10, fc max = 6). The results are 
averaged over 100 independent realizations. 



directed BA networks CBR method can prominently en- 
hance the network synchronizability than the other two 
methods. The suddenly jump of the indicator R for RR 
and DBR methods after adding one link is coursed by 
the existent of a long loop that contains the root node. 
This loop breaks the leadership of root node in driving 
the synchronizing process. 



B. Removing 

When removing links, the CBR method can enhance 
the synchronizability more effectively than the other two 
methods as shown in Fig. [3] In directed regular networks 
[see Fig.[3](a)] , both the RR method and the DBR method 
weaken the synchronizability while the CBR method im- 
proves the synchronizability. By removing links, all three 
methods will destroy the uniform in-degree distribution 
in directed regular networks. However, CBR method can 
generate lots of receptors (i.e., the nodes without any out- 
going link) who are only affected by other nodes. This 
structure is favorable for synchronization since these re- 
ceptors do not send back interruptive information to their 
upstream nodes. When removing the first link from the 
regular network, since the out-degree and the central- 
ity score for all the nodes are the same, it is equivalent 
to a random selection of one link. After that, the node 
whose out link was deleted will have the lowest centrality 
score. Therefore, in the next step a link starting from this 
node will again be removed, until all its out-links were re- 
moved. As a result, this node becomes a receptor. Figure 
@Jb) shows the number of generated receptors of the ob- 
tained networks when cutting 50 links by three methods, 
respectively. The parameter p — corresponds to regular 
networks, p = 1 to random networks, and < p < 1 to 
WS networks. Clearly, CBR method can generate more 
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FIG. 4: (Color online) (a) A simple example to illustrate 
why CBR method has the advantage in shortening the depth 
after removing links from directed BA networks. The nodes' 
centrality scores of the original network are labeled beside the 
nodes, (b) The number of generated receptors when cutting 
50 links from directed WS networks by three methods. The 
original networks are given N = 100 and k — 3 but with 
different p. (c) The depth of the obtained trees by three 
methods from directed BA networks. The original networks 
are given N = 100 and iVo = 10 but with different fc max . 
The results in (b) and (c) are averaged over 100 independent 
realizations. 



receptors than RR and DBR methods. 

In directed WS networks, the CBR method also per- 
forms the best as shown in Fig- EJ^b). At the beginning, 
three methods can enhance the synchronizability by re- 
ducing the maximum in-degree in the networks. After 
the in-degree distribution becomes homogeneous, the sit- 
uation will become similar to the directed regular net- 
works. If we go on removing links by RR method or 
DBR method, the synchronizability will be weaken since 
the homogeneous in-degree distributions are broken. By 
creating many receptors, CBR method can further en- 
hance the synchronizability even the homogeneous in- 
degree distribution gets destroyed. 

In directed BA models, the synchronizability is iden- 
tical in three methods since in acyclic networks the syn- 
chronizability is totally determined by the maximum and 
minimum in-degree. In Fig. [3fc), each stair is corre- 
sponding to a typical maximum in-degree of the net- 
works. Even though the three methods obtain the same 
synchronizability, the convergence time to synchroniza- 
tion are different. By removing the links from nodes in 
relatively lower layers, the obtained networks will have 
shorter depth, and thus converge faster [21J, l22j . Fig- 
ure |U[a) gives a typical example to compare DBR and 
CBR methods. We can see that after removing two links 
by the CBR method from the middle plot the depth of 
the obtained tree is 2. While the DBR method can lead 
to four different trees with equal probability. Since three 



of them are of depth 3 and the rest is of depth 2, the ex- 
pected depth of the obtained tree by DBR is 2.75. The 
statistical results on the depth of the obtained trees after 
removing links from the directed BA models are shown in 
Fig. HJc). Clearly, the CBR method can generate a tree 
with shorter depth than the RR and DBR methods, and 
thus the networks obtained by the CBR method can con- 
verge faster to the synchronized state. See the Kuromoto 
model in Sec. Ill D for numerical evidence. 



C. Rewiring 

An integrated way is to enhance the synchronizability 
by rewiring links while keeping the total number of links 
unchanged. In each step, we firstly remove a link from 
the network according to the strategy of removing links 
(see Sec. Ill B) and then add a link according to the 
strategy of adding links (see Sec. Ill A). Starting from 
different network models, we rewire 300 links via the RR, 
DBR and CBR methods, respectively. The synchroniz- 
ability R of the obtained networks are shown in Table 
HJ The corresponding optimal synchronizability for such 
directed networks is obtained by an evolutionary opti- 
mization algorithm [TH, [l6| . We find that the synchro- 
nizability corresponding to the CBR method is very close 
to the optimal case. In addition, the topology of the ob- 
tained CBR network is very similar to the optimal case, 
such as the hierarchical structures and the existence of 
cores. These results indicate that the CBR method is 
an effective method to enhance the synchronizability via 
link rewiring. 



TABLE I: The synchronizability of the networks obtained by 
rewiring 300 links through RR, DBR and CBR methods. The 
network parameters are set as follows: N = 100 and k — 3 
for regular networks; N = 100, k — 3 and p — 0.1 for WS 
networks; N = 100, k = 3 and p — 1 for Erdos-Renyi (ER) 
networks [JjJ; N = 100, No = 10, fc max = 6 and k = 3 for 
BA networks. The sychronizability of the original networks 
(-Ro) and the result obtained by the evolutionary optimization 
algorithm (i? op t) for the network with N = 100 and k = 3 [l5l ] 
are shown for comparison. 



network 


Regular 


WS 


ER 


BA 


Ro 


156.43 


18.384 


239.49 


6.0000 


Rrr 


3.6975 


3.4459 


3.4043 


3.2231 


Rbbr 


1.3813 


1.4511 


1.4415 


1.3822 


RcBR 


1.2191 


1.2694 


1.3401 


1.3333 


Ropt 


1.17 


1.17 


1.17 


1.17 
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FIG. 5: (Color online) The order parameter of Kuramoto 
model on (a) obtained networks by adding 50 links to di- 
rected regular networks, (b) obtained networks by removing 
50 links from directed regular networks, (c) obtained networks 
by adding 50 links to directed WS networks, (d) obtained net- 
works by removing 50 links from directed WS networks, (e) 
obtained networks by adding 50 links to directed BA net- 
works, (f ) obtained trees by removing links from directed BA 
networks. The network parameters are the same as in Fig. [T] 
The results are averaged over 100 implements. 



D. Kuramoto model 

For the purpose of monitoring the effects of different 
rewiring methods on the convergence time to synchro- 
nization, we study the Kuramoto model 29] in the re- 
sulted networks. The dynamics of N coupled oscillators 
is described by the equation 



K N 



(2) 



where the modulus < r(t) < 1 measures the phase co- 
herence of the population and <p(t) is the average phase. 
r(t) ~ 1 and r(t) ~ describe the limits in which all 
oscillators are respectively phase locked and moving in- 
coherently 

In Fig.[5j we compare the behavior of Kuramoto model 
on the obtained networks from (a) adding 50 links to 
directed regular networks, (b) removing 50 links from 
directed regular networks, (c) adding 50 links to directed 
WS networks, (d) removing 50 links from directed WS 
model, (e) adding 50 links to directed BA networks and 
(f) obtained tree by removing links from directed BA 
networks. 

Generally speaking, the CBR and DBR methods out- 
perform the RR method in the convergence time to the 
completely synchronized state. As mentioned above, the 
emergence of the network core can greatly shorten the 
convergence time. The smaller the core-distance, the 
quicker the whole network converges to synchronization, 
since the core actually drives the whole network to the 
synchronized state. In Fig. [5ja), (c) and (e), it is obvi- 
ous that the CBR and DBR networks have shorter con- 
vergence time than the RR method due to the smaller 
core-distance. 

In Fig. [SJb), the network cannot reach the complete 
synchronization. However, it can be seen that the or- 
der parameter of CBR networks is higher than the RR 
and DBR networks, in accordance with the higher syn- 
chronizability predicted by the master stability analysis. 
Figure [5jd) shows that the networks from the CBR and 
DBR methods converge faster than the ones from the 
RR method. It is because RR is more likely to remove 
long range links, resulting in an increasing of the average 
distance of the network, which usually corresponds to a 
longer convergence time. In contrast, in the directed WS 
networks, the starts of the long range links usually have 
relatively large out-degrees and high PageRank scores, 
and therefore are not likely to be removed according 
to the CBR or DBR methods. As shown in Fig. [SJf), 
the CBR method leads to the smallest convergence time 
among the three methods. As the PageRank centrality 
score can help the CBR method to recognize the hier- 
archical structure of the directed BA model, the final 
spanning trees of the CBR method are always with fewer 
depth than the DBR and RR methods [see also Fig.^c)] 
resulting in the faster convergence. 



where u>i and 9i are the natural frequency and the phrase 
of oscillator i, respectively. The coupling strength K is 
set to be a positive constant, here we fix K — 10 for 
convenience. The collective dynamics of the whole pop- 
ulation is measured by the macroscopic complex order 
parameter 



(3) 



IV. CONCLUSION 

In summary, we proposed a centrality-based recon- 
structing (CBR) method to facilitate synchronization in 
directed networks. The centrality is measured by the 
PageRank algorithm. Numerical analysis on three net- 
work models shows that the CBR method is more effec- 
tive in enhancing the network synchronizability than the 
degree-based reconstructing method and random recon- 
structing method. Specifically, when adding links, the 
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CBR method forms a wisely located core to drive the 
whole network to the synchronized state. When remov- 
ing links, the CBR method generates lots of receptors 
which only receive information and do not disturb the 
upstream nodes. Furthermore, the CBR method can 
also be extended to deal with the link rewiring prob- 
lem. Significantly, the obtained network through CBR 
link rewiring method has very close synchronizability to 
the corresponding optimal case obtained by the evolu- 
tionary optimization algorithm. Finally, the result from 
the Kuramoto model shows that the CBR method has 
advantage in shortening the convergence time to synchro- 
nization. 

In practice, as real systems are frequently manipulated 
by growing new connections and eliminating redundant 
connections [l8j], the CBR method can provide a use- 



ful way to improve synchronizability. In some specific 
systems like ecosystems [HI and neuron systems [43j], 
synchronization should be inhibited. In such situation, 
the inverted CBR method can be applied to decrease the 
synchronizability, namely to add links starting from the 
nodes with low centrality scores to the nodes with high 
in-degrees and to remove the links whose starts have high 
centrality scores while whose ends own low in-degrees. 
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